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ABSTRACT 

Recent surveys have detected Lya emission from z = 4.5 — 6.5 at luminosities as low 
as 10^^ erg s~^. There is good evidence that low numbers of AGN are among observed 
faint Lya emitters. Combining these observations with an empirical relation between 
the intrinsic Lya and B-band luminosities of AGN, we obtain an upper limit on the 
number density of AGN with absolute magnitudes Mb G [—16, —19] at z = 4.5 — 6.5. 
These AGN are up to two orders of magnitude fainter than those discovered in the 
Chandra Deep Field, resulting in the faintest observational constraints to date at 
these redshifts. At z = 4.5, the powerlaw slope of the very faint end of the luminosity 
function of AGN is shallower than the slope observed at lower redshifts, /3i < 1.6, at the 
98% confidence level. In fact, we find marginal evidence that the luminosity function 
rises with luminosity, corresponding to a powerlaw slope /?/ < 0, at magnitudes fainter 
than AIb ~ —20 (75% confidence level). These results suggest either that accretion 
onto lower mass black holes is less efficient than onto their more massive counterparts, 
or that the number of black holes powering AGN with Mb > — 20 is lower than 
expected from the A/bh — f relation by one-two orders of magnitude. Extrapolating 
from reverberation-mapping studies suggests that these black holes would have A/bh = 
10^ — 10 ''M0. To facilitate the identification of AGN among observed Lya emitters, 
we derive observational properties of faint AGN in the Lya line, as well as in the X-ray 
and optical bands. 
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1 INTRODUCTION 

The quasar luminosity function (QLF) describes the space 
density of Active Galactic Nuclei (AGN) as a function of 
luminosity and redshift. The QLF encodes information on 
quantities like the black hole number density per unit mass, 
and the gas accretion efficiency. It therefore constrains phys- 
ical models of AGN and of super massive black hole forma- 
tion. The optical (B-band) QLF at z ^4 has been deter- 
mined accurately for luminosities^ exceeding Lb J>1O^^I/s,0 
from the 2dF quasar survey (which corresponds to abso- 
lute magnitudes Mb < —22, Boyle et al. 2000; Groom 
et al. 2004). At higher redshifts the optical luminosity 
function of luminous quasars has been determined from 
quasars in the Sloan Digital Sky Survey jFan et al.ll20ol 
at Mb ~ —27. In addition, deep Chandra and XMM imag- 
ing has constrained the X-ray QLF at X- Ray luminosi- 
ties as low as Ljs 
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the X-Ray QLF to constrain the B-band QLF down to 
fainter optical luminosities {AIb ~ —22) at z ^4. At these 
redshifts, no observational constraints exist on the optical 
QLF at fainter luminosities. However, the details of the ex- 
istence and evolution of low luminosity AGN is crucial for 
our understanding of the growth of low mass black holes, 
and of their role in the formation of super massive black 
holes. 

In this paper we demonstrate how existing Lya surveys 
may be used to constrain the B-Band QLF at absolute mag- 
nitudes as low as Mb = —16. Existi ng wide-field narrow- 
band surveys (e.g. iRhoads et allbOOOD are optimised to de- 
tect Lyof line emission from high-redshift galaxies, and in 
order to maximise their detection r ate, deeply image field s 
as large as 1 deg^ on the sky (e.g. iTaniguchi et aljl2005t) . 
The combination of wide field and deep images, together 
with concentration on a strong emission line, rather than 
continuum, allows these wide-field narrow-band surveys to 
put stringent constraints on the number density of AGN 
with Mb ~ —19. Constraints on even fainter AGN are de- 
rived from deep spectroscopic surveys of regions around in- 
termediate redshift clusters of galaxies, which offer an ultra- 
deep view into the high redshift universe (z=4.5-6.7) through 
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strong gravitational lensing (with magnification factors of 
10-1000, Santos et al, 2004). 

Throughout this paper, the terms AGN and quasar are 
interchangeable and refer to broad lined active galactic nu- 
clei, also known as 'Type F AGN. The outline of this pa- 
per is as follows: In §|21we relate the Lya luminosity of a 
quasar to its B-band luminosity. In § |3 we summarise con- 
straints on the number density of Lya emitters at high red- 
shift, and discuss the abundance of quasars among these 
sources. We show how this existing data constrains the faint 
end of the quasar luminosity function. In §2] we calculate 
the observable Lya properties of AGN and use these to 
physically interpret our results. We also model the appear- 
ance of faint AGN in the optical and X-ray bands. Finally, 
in § |5] we discuss the possible cosmological implications of 
our work, before presenting our conclusions in § |S| We use 
the VKM^P cosmological parameters: {i}m,^A,^b,h,YH_s) 
=(0.3,0.7,0.044,0.7,0.24) llSpergeret al.1 (200311 throughout 
the paper. 



2 THE RELATION BETWEEN LYa AND 
B-BAND LUMINOSITIES IN AGN. 

There is a very simple relation between the Lya and B- 
band luminosities of a quasar. The continuum flux density 
of a quasar may be denoted either by Fe ( in units of erg 
cm~^ erg ~^) or Fx ( in units of erg s~^ cm~'^ A"'^)- 
These quantities are related by the equality XFx = EFe, 
where E — hpc/X is the photon energy and hp is Planck's 
constant. Assuming the Lya line has an equivalent width 
EW, the Lya luminosity of a quasar becomes 

/ EW \ 

Uyc = EWFx\ ° = ( ,-, ][\Fx] (1) 

^ '''a=1200A V 1200 A/ ^A = 1200A ^ ' 

The B-band luminosity can be estimated from 



Lb = BW Fx 



(S)t^-^^^A=4400i' (2) 



where we assumed the B-band filter to be centered on A = 
4400 A with a width of BW= 800 A. The ratio of L^ya and 
Lb becomes 
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have computed the characteristic 



ISazonov et"ai] ll200. 

angle- averaged, broad-band spectral energy distribution of 
the typical quasar. They show that in the range E = 1 — 10 
eV, EFe oc E'^ , with 7 = 0.4. The exact va lue of 7 varies 
between individual objects. iFan et al] ll200 found that at 
z ~ 4.5 the spread in 7 could be represented by a Gaussian 
with a-y ~ 0.3 in their sample of 39 quasars. 

The intrinsic Lya equi valent width for AGN is uncer- 
tain. ICharlot fc FalJ Jl993l) showed that AGN which are 
completely surrounded by neutral hydrogen gas have Lya 
EWs of 827a-i(3/4)" A, where a is the spectral index 
blueward of the Ly a line. According to the template of 
ISazonov et"ai] tOOi) . a = 1.7, which yields EW ~ 300 A. 
However, the observed Lya EWs of bright AGN are typically 
in the range 50-150 A, which co uld be attributed to dust 
attenuation of the Lya line (see ICharlot fc Falllll99j . and 



references therein) and/or scattering in the IGM (see 14.111 . 
iFan et all ll200l[) represent the distribution of observed EW 
by a Gaussian centered on EW= 50 A, with (Jew = 14 A. 
^ In § 14.11 and Appendix ^ we have argued in detail that 
~ 50% of all Lya photons emitted by AGN is transmit- 
ted through the IGM (at z < 6; at z = 6.5 this number is 
^ 30—40%). If half of the emitted Lya photons are scattered 
in the IGM, then the intrinsic distribution of EW should be 
centered on EW= 100 Awith ctew = 30 A. The distribution 
of the ratio L^ya/Ls was obtained via Monte-Carlo, under 
the assumption that: 1) the slope 7 follows a Gaussian distri- 
bution with 7 = 0.4 ± 0.3; and 2) the Lya equivalent width 
follows a Gaussian distribution with EW=100±30 A, which 
is truncated at EW= and 600 A. The ratio Li^ya/ Lb and 
it's 95% confidence interval are found to be 
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Thus, within the uncertainties we may approximate the Lya 
luminosity of an AGN as being equal to its B-band luminos- 
ity. 

It is important to stress that the results that are de- 
rived in this paper assume that this ratio does not change 
toward lower AGN luminosity. This eflFectively means that 
we assume the AGN's Lya EW is independent of luminos- 
ity. However, observations of both high and low redshift 
quasars and Seyfert I galaxies revealed that the EW of the 
CIV 1549 line increased toward lower luminosities. The in- 
crease of the EW in certain emission lines in spectra of AGN 
toward lower AGN luminosities is known as the 'Baldwin- 
Effect'. T his effect, however, has not b een observed for the 
Lya line jMushotzkv fc Ferlandl ll984>l. We point out that 
even if the Baldwin-Effect were applicable, then this would 
increase the ratio LLya/is. For a fixed Lya luminosity, we 
would then be able to probe AGN with fainter Lb, which 
would strengthen the main results presented in this paper. 



3 CONSTRAINING THE FAINT END OF QLF 
USING LYa SURVEYS. 

3.1 The Number Density of AGN at z = 4.5 - 6.5. 

Several Lya surveys have derived the number density, 
riLya, of Lya emitters brighter than a minimum de- 
tectable Lya luminosity, Llv^.c. Examples includ e surveys 
at z = 4.5 iD awson ct al. 2004*) , z = 5.7 (e.g . iHu et alJ 
I 2OO4I: lOuchi et al 2005.: .Shimasaku et al.l 1200611 . 2 = 6.5 

I Taniguchi et al J 120051: iKashikawa et alJbOO^T and z > 4.5 
Santos et all2004^ . Results from surveys that were used for 
this paper are summarised in columns 1-4 in Tabled 

Since the detected Lya emitters are generally identi- 
fied as galaxies, the number densi ty of quasars w i th Ll vq > 
L-Lya.c must be smaller than WT,,,fv. Dawson et al.l (I2OO4II find 
that no AGN are among the Lya emitters in their sam- 
ple in a total comoving volume of 1.5 x 10® Mpc'^. This 
finding was based on the absence of broad line emitters 



^ iFan et all i200 j) actually measure the sum of the Lya and NV 
equivalent width to be EW= 70 it 20 A. The Lyo line makes up 
~ 70% of the total EW. 
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Table 1. Lya Survey's Constraints on nLya and the corresponding upper limits on 'I'(I/B,mini z)- 



Z /min 


-^LycKjC 


n-Lya "AGN 


^B.min 


^(iB.min.z) 


(10~^* erg s~^ cm~2) 


(10*2 erg s-l) 


(10-* Mpc-3) (lO"'' Mpc-3) 


10SLb,q 


(lO-f' Gpc-3 L-1) 


4.5^ 16 


6 


1.6 < 2 


20 


< 0.08 {/3 = 1.1) 


5.72 6 


4 


3.9 < 58 


14 


< 2.6 (/3 = 1.6) 


6.6^ 4.1 


6 


1.2 < 35 


20 


< 1.1 (p = 1.6) 


> 4.5"* 0.1 


0.06 


102.0 < io*-8 


0.7 


< 104-9 (/3 = 1.6) 


> 4.5* 0.3 


0.2 




2.2 


< lO^ '' (/3 = 1.6) 


(1) Data fromlDawson et al.l J2004 


\j. Note that the actual detection limit is quoted as r 


- 2 X 10-1'' erg s" 


1 cm— 2 We have slightly 



revised their detection limit downwards, as 27% of their spectroscopically confirmed LAEs have a flux < 2.0 X 10-^' erg s— ^ cm~^ , 
with a mean flux of ~ 1.1 X erg s-^ cm-^. As argued in S 13.11 the 5 — cr detection threshold for AGN lies higher 

by a factor of 1.4 at ~ 1.6 X 10— erg s-^ cm-^. 

(2) Data froml Ouchi et al.l ^2005^ 

(3) Data from lTaniguchi et al.l ^2005^ and iKashikawa et alJ i2006^ 

(4) From the deep spectroscopic survey of gravitationally lensed regions bv lSantos et al.l j20o3). We have increased Lb, min by a 
factor of vTO (relative to that given by Eq|l} to account for the reduced sensitivity to detect the AGN's broad Lya emission 
lines in this spectroscopic survey (see the main text for more discussion on this). 



in both their narrow band survey and in follow-up high 
resolution spectra of 18 confirmed z — 4.5 Lya emitters. 
Furthermore, these objects lacked high-ionization state UV 
emission lines, symptomatic of AGN activity. Additionally, 
deep X-Ray observations of 101 Lya emitters by Wang et 
al. (2004, also see Malhotra et al. 2003), revealed no X-ray 
emission from any individual source at a 3 — ct detection limit 
of Jb,5-io.okeV = 6.6 X 10~^^ erg s~-^ cm~^ (but this may 
not be that surprising as we discuss in § 14.31 . The absence 
of AGN in the survey volume yields an upper limit on the 
number of AGN, A'^agn ^ 3 at the 95% confidence level.* 

Similarly, Quchi et al. (2005) found 515 objects within 
their survey volume of 1.3 x 10^ Mpc"^ that are probable 
z = 5.7 Lya emitters. Follow-up spectroscopy of 19 emitters 
revealed that their spectra are too narrow (wfwhm <,500 km 
s-^) to be AGN. If the sample of 515 detections contains a 
fraction of AGN, /agn ;^0.15, then the probability that a 
random sample of 19 contains no AGN is (1 — /agn)^^ Jj0.05. 
This yields a 2 — ct upper limit on the number of AGN, 

Na q-n <;515fAGN =7 5- 

iTaniguchi eHU J2005ft and IKashikawa et al.l ll200ell 

The survey volume of lDawson et al] j20o3) is determined by 5 
overlapping filters each of width 80 A, yielding a total wavelength 
coverage of 240 A. As will be demonstrated in 5 14.11 the observed 
FWHM of the Lya line of AGN at this redshift is typically ~ 60- 
70 A, which implies that AGN would have been de tected in more 
than two filters. From narrow- band imaging alone. IPawson et alj 
i2004h find no evidence for the presence of broad line emitters. 
If we represent the filter transmission curves with top hats of 
width 8OA, and an AGN's Lyo emission line with a Gaussian 
with a FWHM of 65 A, then on average ~ 70%, 60% and 40% 
of the AGN's total Lyo flux is detected by three adjacent filters. 
This implies that 1) for an AGN to be detected in 1 filter at the 
5-cr detection threshold, requires its Lya flux to be 1/0. 7~ 1.4 
times higher than the surveys detection threshold, and 2) for an 
AGN to be detected in 3 adjacent filters requires its peak fiux to 
lie within the inner three narrow band filters. This reduces the 
survey volume in which these AGN can be identified by a factor 
of (240 - 2 X 40)/240 = 2/3. 

* This is because for a Poisson distribution P(n,,/i) = e-''/i"'/n! 
with an expected number of events p. = 3.0, the probability of 
having more than events occur is 95%. 



found 53 candidate z = 6.5 Lya emitters within their sur- 
vey volume of 2 X 10^ Mpc'^ centered on the Subaru Deep 
Field. Follow-up spectroscopy of 22 candidates revealed a 
combined 17 objects that are likely to be z = 6.5 Lya emit- 
ters. The line widths of their objects lie in the range 180-480 
km s-i. This, in combination with the lack of NV A = 1640 
A emission, strongly suggests no AGN are in their sample. 
Following the discussion above, if the sample of 53 candi- 
date Lya emitters contains a fraction of AGN, /agn ^0.12, 
then the probability that a random sample of 22 contains 
no AGN is (1 - /agn)^^ ^0.05. This yields a 2 - a upper 
limit on the number of AGN, A'agn ;^53/agn ~ 7. 

In their deep, blind, spectroscopic survey of regions that 
utilized strong-lousing magnification by 2 = 0.2 clusters of 
galaxies, Santos ct al. (2004) found 3-5 Lya emitting objects 
at 2 > 4.5 at unlensed flux levels as low as ~ 3 x 10-^^ erg 

cm- 2. The l ine wi dths of these sources are not reported, 
but lEUis et alJ i200ll) found the line width of one z = 5.6 
Lya emitter to be ;^100 km s-^, and is therefore unlikely 
an AGN. However, to be conservative, our upper limits are 
derived under the assumption that all these Lya emitters 
could be AGN. 

For each survey, we obtain upper limits on the num- 
ber density of AGN from tiagn ^A^'AGNVgurvcy The detec- 
tion limit of each survey, /min, is converted to its mini- 
mum Lya detectable Lya luminosity, LLya,c, using L^ycc — 
4-Kdj;^{z) /iim/(T). Here, dL{z) is the luminosity distance to 
redshift z and T is the mean transmitted Lya flux from the 
AGN through the IGM. As will be discussed in § 14.11 we as- 
sume that for AGN T = 0.5 when 2 < 6 and T — 0.3 when 
z — 6.5. Valu es of /lim are given i n Table Q For the survey 
performed bv lSantos et al.l ^2004^ ■ the calculation of uagn 
is far more complicated, since both the detection limit and 
survey volume are functi ons of the po s ition on the sky and 
of redshift. Fortunately, ISantos et alJ (|20o3) fully account 
for this and provide number densities of Lya emitters with 
L > Z/crit, for several values of Lcrit- We obtain upper limits 
on riQso from the upper boundary of the 95% confidence 
levels on their quoted nLya(i > Lait) for log[I/crit] =40.5 
and 41 (their Figure 12. Our LLya.c = 2Lcrit because we 
account for 50% loss of flux in the IGM). Because the Lya 
emission lines of AGN are expected to be much broader 
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than those of galaxies (typically by a factor of ~ 10, see 
§ 14. in . these are more difficult to detect in spectroscopic 
surveys, in which the AGN's Lya flux would be spread out 
over ~ 10 times as many frequency bins. This would in- 
crease the total noise by a factor of ~ y^lO). To represent 
the decreased sensitivity to Lya emission lines emitted by 
AGN, the value of LB.min shown in Tabled was obtained 

from LB.min = V^10)LLya,c/0.7. 

A summary of the minimum detectable Lya luminosity, 
L-Lya,c and our constraints on the number densities of AGN 
at z — 4.5, 5.7 and 6.5 is given in Tabled By converting 
L-Lya,c to a minimum B-band luminosity I/B,miii (Eq. EJ, 
we are able to constrain the B-band QLF. This process is 
described in the following subsections. 



3.2 Using the Number Density of AGN at 
z — 4.5 — 6.5 to Constrain the QLF. 

The luminosity function '^'{Lb, z)dLB is defined as the num- 
ber of quasars per unit comoving volume having rest-frame 
B-band luminosities in the range between Lb and Lb +d,LB 
at redshift z. The number density of AGN with Lb > LB.min 
is then given by 



nAGN{LB > Lb 



^{LB,z)dLB (5) 



To constrain the luminosity function at I/B,miii we must 
know the functional form of ^'(Ls,^), which is not con- 
strained observationally at the redshift range of interest, 
z G [4.5, 6.5]. We first constrain '^'[Lb, z) at I/B,min at z > 5 
under the assumption that the faint end of the quasar lu- 
minosity function observed at lower redshifts applies here 
as well. This is followed by a self-consistent constraint on 
V['(Ls, z) and its slope at z — 4.5. 

3.2.1 Constraints on \['(L_b,z) at z > 5. 

The following double power law provides a good represen- 
tation of the observed quasar luminosity function at z ^ 3 
jBovle et al.lE)0al 



*(Ls,2) 



*./L*(2) 



( Lb \l3h , ( Lb \^ 
\L*(z)) \L*(z)) 



(6) 



The slope at the bright and faint end of the luminosity func- 
tion is Ph = 3.52 and Pi = 1.66, respectively. All redshift 
dependence lies in the transition luminosity L,{z). If the 
same parameterization holds at z > 5 and I/B,min ^ Lt{z), 
then the integral in Eq (|SJ is dominated by the faint end of 
the luminosity function and the /3h— term may be omitted. 
Eq ^ then becomes 



nAGN(I/fl > Lb.u 



^B.min 



Y*(iB,min, z) 



(7) 



which shows that the upper limit on nAGN(is > ifi.min; z) 
translates directly to an upper limit on ^'(I/B,min, z). These 
constraints are summarised in the last column of Table 
(at z = 5.7 and z — 6.5). In Figure Q we compare the above 
constraint on the quasar B-band luminosity function (de- 
noted by the filled circles at logfLs/I/B,©] ^10 in central and 
right panel), with data at higher luminosities. Each panel 
corresponds to the redshift printed in the upper right cor- 
ner of the figure. The open triangles at log[Lfl/LB,©] > 13 



represent data from Fan et al (2001b). The open squares 
are de rived from the X-Ray data presented bv lBarger et alJ 
i2003l) . who plot the number density of AGN as a func- 
tion of redshift in the rest frame soft X-Ray (E=0. 5-2.0 
keV) luminosity range Lx = 10*'' — 10*'' ergs s~^ (here- 
after, the 'faint bin') and Lx = 10** - lO*'^ ergs s~^ (here- 
after, the 'bright bin') . For the ISazonov et al.l i2004l) tem- 
plate spectrum, the ratio of the rest frame soft X-Ray (de- 
noted by Lx) to blue band luminosity is Lx / L b ~ 0.5. 
The faint and bright bins of lBarger et alJ ll2003h therefore 
constitute the range Lb = 5 x 10^" — 5 x IO'^'^Ls^q and 
Lb = 5 X lO'^'^ — 5 X W^^Lb.q, respectively. To convert the 
X-Ray number density to the luminosity function, '^'{Lb, z), 
we simply divided rifaint by the total luminosity width of the 
bin. We use nfeint = i^-^'ltl) ^ 10 Mpc"^ in both bins. 

The error-bars on our constraints reflect the uncertainty 
in the exact ratio of Lb / iLya (Eq. 2J and denote the 95% 
confidence levels. Note that upper limits on >I'(LB,min, z) 
as a function of LB.min have power law slopes of —1 (see 
Eq. |7|l. The constraints from the wide fiel d surveys lie on 
or slightly below the model predictions of IWvithe fc Loebl 
(2003) (shown as the solid lines). This implies that the ab- 
sence of z = 5.7 and z = 6.5 AGN in wide-field narrow-band 
surveys rules out this model at > 95% confidence levels at 
both redshifts, at these low luminosities. 



3.2.2 Constraints on'^{LB, 



at z ^ 4.5. 



The above constraints at z > 5 assumed that /3; = 1.66. 
However, as will be discussed below, at z = 4.5, extrap- 
olation of this powerlaw from the X-Ray constraint at 
log[L_g/L_g,Q ~ 11] to I/B,min, would havc resulted in a num- 
ber density of AGN that lies above the implied upper limit. 
To better constrain the luminosity function, ^'(LB.min, z), 
and its slope at z = 4.5, we combine the upper limit on the 
AGN number density from the Lya survey w ith the AGN 
number densities derived from X-ray data jBarger et alJ 
!2003; Co wic et al. 2003) . We use nbright = (1.0±0.5) x 10"® 
Mpc"^ and nfai„t = (3.3 ± 1.3) x 10"*^ Mpc"^, which we 
obtained by i nterpolating between the z = 3.5 and z = 5.7 
data points of lBarge~ et al.U2003l) . 

We assume the luminosity function, ^'(L_b,z) to be 
of the form kL~^^ . The number density of AGN in the 
faint X-ray bin, nfaint, fixes fc as a function of Extrap- 
olating this luminosity function to lower luminosities, gives 
us the expected number density of AGN in the Lya bin, 
magn = k f^^'""^'' dl, L~^' . One of the selection criteria 
iDawson et alJ (|200^ use to select z = 4.5 Lya emitter can- 
didates states that a z = 4.5 candidate must not be de- 
tected in the Bw band of the NOAO deep wide-field survey 
IMalhotr a & Rhoads 200^). We show in § |01that this im- 
plies that Lb, max ~ 2 x 10^^Lb,q. For a fixed we obtain 
a unique tiagn. Since tiagn is only known as an upper limit, 
this allows us to put an upper limit on /3; as follows. 

In the absence of AGN within the volume probed by 
the Lya survey, the probability that the true number lies 
between n and /j. + d/j, is given by dP =exp(— ^)d/^ (where 
we have assumed a Poisson distribution) . This may be recast 



P(< m) = 1 - e" 



(8) 
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Figure 1. Ly« constraints on the quasar B-band luminosity function (QBLF), >I'{Lfl, z) at redsliifts z = 4.5 (ie/t panel), z = 5.7 {middle 
panel) and 6.5 {right panel), plotted as filled circles at logfLs/iB,©] ^10. Tlie error-bars on the Ly« constraints denote 95% confidence 
limits. The open circle at z = 4.5 denotes the upper limit o btained from spectro scopic data alone 13.31 . The SDSS and X-Ray data 
are taken from Fan et al (2001,2001b) {open triangles) and iBarger et all i20n3l) {open squares), respectively. This figure shows that 
existing Lya surveys allow the QBLF to be constrained at B-band luminosities as low as logfLs/LB,©] = 8.5 (or at absolute magnitude 
Mb ~ —16). Our constraint at z = 4.5 suggests that the QBLF fiattens for logfLs/Ls,©] <^11 (also see Fig. 1^. 



which gives the probability that the expected number of 
AGN within the survey volume is less than ju. The expected 
number of AGN is given by 

^^ = Vsurvoy / dL ^(L.z). (9) 

For a fixed Lb, max, we find that /i is a function of LB,min, 
?^faint and /3; only. In the left panel of Figure |21 the thin 
black solid line shows the probability P(< Pi) as a func- 
tion of Pi, for the fiducial values of I/B.min = L/Lya.c and 
rifaint = 3.3 X 10"^ Mpc-^ In this case P(< 0.72) = 0.95, 
i.e. the slope is /3; ;^0.72 at the 95% confidence level. To il- 
lustrate the dependence of P{< f3i) on LB.min, and nfaint, we 
have also plotted P{< Pi) for (nfaint, Lb, min) = (2.0 x 10~® 
Mpc"^, LLya,c) (red dotted line) and (3.3 x 10"*^ Mpc"^, 
0.5LLya,c) {blue dashed line). The 2 — a upper limit on /3; 
increases with decreasing rifaint and increasing I/B,inin. To 
eliminate the dependence of P(< Pi) on I/B.min and rifaint, 
we marginalise over these parameters 

-r, , f , f ,r dP dP 

■Pmarg(< Pi) = arlfaint / dLB.min-: -JZ X 

J J "rifaint ai>B,min 

xfl-e""*"'-""-^'^'""-''')], (10) 

where dP/drifaint and dP/dLB,min are the probability dis- 
tributions for rifaint and LB,min, respectively. We choose 
dP/drifaint to be Gaussian with rifaint = (3.3 ± 1.3) x 
lO"'' Mpc~^, and dP/dLB,min to be log normal with 
log(I/B,min/Ls,0) = 10.3 ± 0.2 (this range is motivated by 
the distribution in the ratio LLya/Ls, found from Eq. |l|in 
§|2|l. The marginalised constraint is given by the thick solid 
hue in Figure |5| which shows that /3; ;^1.1 at the 95% confi- 
dence level. The right panel of Figure |5| shows the expected 
value of *I'(I/B,min, a) as a function of Pi. The constraint 
Pi < 1.1 yields log[\l/(I/B,min, z)] < —7.1 and is shown as the 
filled circle. This 2 — cr upper limit on ^'(I/B.min, 2) at 2 = 4.5 



is also shown in the left panel of Figure Q as the filled cir- 
cle. The open triangles at \og[LB / Lb,q] > 13 represent data 
from Fan et al (2001) and the open squares represent the X- 
ray constraints derived from the data presented by Barger 
et al. (2003) and Cowie et al. (2003). Figure [T] implies that 
at 2 = 4.5, the QLF flattens below \og{LB / Lb,q) ~ H, as 
may also be seen via the model quasar luminosity function 
from Wyithe & Loeb 2002 and 2003 (plotted as the black 
solid line in Fig. 0. It is worth pointing out how the ex- 
act marginalised constraint on Pi depends on the assumed 
survey sensitivity. If we decrease/increase /um by a factor 
of 2 then /3; < 0.8/1.4 at the 95% confidence level, respec- 
tively. The value of ■Pniarg(< 0.0) barely depends on /min. 
The decreased sensitivity corresponds approximately to the 
sensitivity limit reported by'-Dawson ot al. (2004), while the 
increased sensitivity corresponds to the sensitivity derived 
from their actual faintest Lya detections. 

The constraints obtained so far only used the non- 
detection of broad line emitters in the narrow band survey. 
Dawson et al. (2004) also obtained spectra of 25 of their 
candidate LAEs. Below, we constrain 'i'{LB,z) using this 
spectroscopic data. 

3.3 Constraints on '^{Lb,z) at 2 = 4.5 from 
Spectroscopic Follow-Up. 

iDawson et alJ i2004|) obtained follow-up high resolution 
spectra for 25 candidate LAEs. Of these 25 candidates, 18 
were genuine 2 — 4.5 objects and the other 7 were either 
not detected (6) or a lower redshift [O II] emitter (1). All 
18 2 = 4.5 objects were identified as galaxies, based on 1) 
the observed narrow physical line widths of the Lya lines 
(At; <;500 km s~^, also see § 14. IL and 2) the lack high- 
ionization state UV emission lines, symptomatic of AGN 
activity, in their spectra. If the sample of 350 detections 
contains less than 10% AGN, /agn ;$0.1, then the prob- 
ability that a random sample of 25 contains no AGN is 
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/Sj=-aiog ^/d\og Lg /3j=-aiog */aiog Lg 

Figure 2. Left Panel: The probability P(< /3;) that — 9 log ^/O log Ls is less than /3j. The i/iin black solid line represents the constraint 
for the fiducial model with Le.min = ^Lya.c nf^^^^ = 3.3 X 10~® Mpc~^, for which P{< 0.72) = 0.95. The effect of lowering Uf^^^^ {red 
dotted line, ^faint — 

2.0X 10"'' Mpc-3) and Lb, min {blue dashed line, LB,min — 0.50ZyLycK^c) is to reduce and increase -P(< respectively. 
The thick solid line shows the probability, Pmarg(< /S;), after marginalizing over Lb, min and rifaint. We find marginal evidence (75% 
confidence level) that the luminosity function falls towards lower luminosities (Pmarg(< 0.0) = 0.75) below log[Ls/Ls q] = 11. Right 
Panel: The expected value of 'I'(LB,mini 2) is plotted as a function of /3; for the models shown in the left panel. The marginalised 2 — a 
upper limit (also shown in Fig0 is shown as the large filled circle, which implies log['I'(LB mini 2)] < —7.1. The light and dark grey 
areas demarquate the region of (/3; , \E')-space excluded at the > 68% (l-O") and > 95% (2-(t) confidence level, after marginalising over 
"faint and Lb, min- For reference, the model of lWvithe' fc Loebl ^2003^ predicts log"!- ~ -5.7. 



~ (1 - fAG^f^ ^0.05 (In §|01we show that /agn ^0.1 
is consistent with the constraints imposed by the X-Ray ob- 
servations of Wang et aL 2004) . From the spectroscopic data 
alone, we can therefore put a 2 — a upper hmit on the num- 
ber of AGN within the original survey volume (1.5 x 10^ 
Mpc^) at Nacn ^350/agn = 35 down to a Lya flux of 
1.1 X 10~^^ ergs s~^ cm~^. Using Eq Q we obtain an up- 
per limit on ^'(Lfl,2), which is shown as the open circle in 
Figure This constraint on 'I'(Z/b,2) is weaker than our 
original constraint, but still rules out the model luminosity 
function at > 95%. 



4 THE NATURE OF VERY FAINT AGN AND 
THEIR OBSERVABLE PROPERTIES. 

We have demonstrated how Lya surveys may be used to 
constrain the (very) faint end of the quasar luminosity func- 
tion. This constraint is purely empirical and we stress that 
so far no details of the process of gas accretion onto black 
holes have been assumed. The following discussion focuses 
on the nature of faint AGN (with log[Ls/LB,0] <10). First, 
we present a more general calculation of the observable Lya 
properties of AGN. The result of this calculation allows us 
to determine the range of black hole masses probed by Lya 
surveys and may be used to identify AGN among known 
Lya emitters, especially in combination with our estimates 
for the observable properties of these faint AGN in the op- 
tical (§|3iU and X-ray (§|3inj bands. 



4.1 Observable Lya Properties of Faint AGN 

First we discus s the intrinsic Lya properties of AGN. The 
Kcispi relation iKaspi et aLlbOOCt iPeterson et al.ll20o3) re- 
lates the mass of the black hole powering the quasar, Mbh, 



to its continuum luminosity at 5100 A 

MBH^7.6±1.3X10^(^^M^)°-Me. 

V 10**erg s^ / 



(11) 



We define the ratio 7^ = [AFa(1200 A)]/[AFa(5100 A)] = 
[^ ^^(1200 A)]/ \EFe(5W0 A)] (see The template given 
by ISazonov et al. (2004) yields TZ = 1.8. For a given Lya 
EW, the Lya luminosity is therefore uniquely determined 
by Mbh through 



T 1 « ^ EW N 



Mbh 



VlOOA/ V7.6 X lO^M© 



(12) 



where LLyQ,,43 is in units of 10'*'' ergs s ^. The Lya emis- 
sion lines from A GN are broader than tho se of galaxies. 
iKaspi et all (|200(t and lPeterson et al.l lj2004h found the Half 
Width at Half Maximum (HWHM) of the Balmer lines to be 
related to the continuum luminosity at 5100 A. If the same 
relation holds for the Lya line, then its HWHM is given by 

Mb 



.HWHM = 2 X 10^ ( ^ ^^-^^J kms-. (13) 

Note that the observations presented in lKaspi et alJ (l200ff) 
(their Figure 7) suggest that vhwhm does not increase be- 
yond uhphm ~ 3 X 10"^ km s"'^ at low masses. 

To obtain the AGNs' observed properties in the Lya line 
from the above relations, we require knowledge of the frac- 
tion of transmitted Lya-flux through the IGM as a function 
of frequency. We first focus on the observed Lya properties 
of AGN a,t z = 4.5 and z = 5.7. Although the universe is 
fully reionised at 2: < 6 (e.g. iFarTe t al. 200^, a trace quan- 
tity of neutral hydrogen is sufficient to scatter Lya blueward 
of the Lya line center out of our line of sight. For this reason 
the IGM, to first order, erases the blue half of the Lya line, 
so that the IGM transmission is T = 0.5. Various refine- 
ments of this scenario are possible. Infall of the IGM around 
massive objects causes the IGM to erase a part of the red 
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side of the Lya line as well iBarkanall2004l) . which reduces 
T. On the other hand, the proximity effect around these 
fainter AGN increases T. These two effects counteract each 
other, and ignoring both does not add a significant error 
to our estimate of T. This is especially true when consider- 
ing the large width of the Lya line emitted by AGN. When 
the IGM erases the blue half of the Lya line, the observed 
fHWHM is reduced by a factor of 2 relative to the value in 
Eq (1131 . Note that dhwhm = 1500 km s~^ corresponds to 
an observed FWHM of (1 + 2)2uhwhm/c ~ 66A at z = 4.5 
(this value was used in S 13.11 . Note that one has to be care- 
ful not to confuse the observed FWHM of the Lya line with 
its EW (it may be particularly confusing since the observed 
rest-frame EW of AGN are 50 A, comparable to the value 
of the observed FWHM here). 

At 2; = 6.5, the IGM is believed to contain a signif- 
icant fraction of neutral hydrogen IjW yithc & Locb 20o3: 
iMesineer fc Haima n'2004': ' Fan et al.l2005() . For Lya sources 
embedded in a neutral IGM, the damping wing of the Gunn- 
Peterson trough can extend to the red side of the line and 
erase a significant fraction of the total Lya flux. However, 
this effect is re duced when the Lya source is surrounded b y 
an HII region JCen fc Haimanll2000l : (Madau fc Reesll200Cl) . 
Moreover, sources with sufficiently broad emission lines 
(I'HWHM > 300 km s~^) can remain detect able even in the 
absence of such an HII region (lHaiman ^2002''). In Appendix|Xl 
we show that a representative number for the IGM trans- 
mission T for faint AGN embedded in a neutral IGM at 
z — 6.5 is T = 0.3 and that the observed uhwhm is reduced, 
again by a factor of ~2 relative to the value in Eq l|13|l . Fur- 
thermore, we show that the Gunn Peterson damping wing 
may cause the observed line center of AGN at 2: = 6.5 to be 
redshifted by up to ~ 2000 km s~^ relative to other emission 
lines (see Appendix EJ- 

Next, we write the total detectable Lya flux of AGN as 

/ T w EW \ /I -Hz\-2-8 _i _2 

in which T ^ 0.5 for z < 6 and T ^ 0.3 for z = 6.5. 
To allow fa to be written as a simple function of {1 + z), 
we approximated the luminosity distance by dL{z) = 3.6 x 
10*([1 -I- z]/5y-* Mpc, which is within < 6% of its actual 
value in the range z — 3 — 10. Eq 1141 shows that existing 
Lya surveys could have detected AGN powered by black 
holes with masses of Mbh ~ 10*^ - lO'^M© at z = 4.5 - 6.5. 

Equation 1141 assumed that the Kaspi-relations derived 
from observations of luminous AGN also relate the black- 
hole mass to the Lya luminosity for low luminosity AGN. 
Alternatively, if we assume that faint AGN are powered by 
less massive black holes accreting at their Eddington limit, 
then these black-holes are up to ~ 5 orders of magnitude less 
massive than those of the observed SDSS quasars, yielding 
Mbh ~ a few 1O*M0. Combined with the above argument, 
we therefore conclude that existing Lya surveys could have 
detected AGN powered by black holes with masses in the 
range Mbh = 10'' - lO^M© at z = 4.5 - 6.5. Here, the lower 
end of this mass range is probed only by the spectroscopic 
surveys of gravitationally lensed regions. 



4.2 AB-magnitudes of Faint AGN 

Using the t e mplat e spectrum for AGN given by 
ISazonov et al.l i2004l) . we calculate the apparent AB- 
magnitude of faint AGN as a function of observed 
wavelength. The AB-magnitude of an object with flux 
density F„ (in ergs s~^ cm~^ Hz^') is 

mAB = -48.6 - 2.5 log(F,). (15) 

At rest frame energies -E < 10 eV and E > 10 eV, the 
AGN continuum follows the power-laws Fi, oc and 
oc i^~''^, respectively JSazonov et alJl2004D . For an AGN 
whose Lya flux is /a, 17 x 10~'^ ergs s~' cm~^, the apparent 
magnitude may then be written as 

mAfl(A) = 26.4-2.5[log(^^) +log/C(A)], 

where 

- (¥)(w)(¥)(xr:)°"-v. 

- <o(^r)(w)(f)(^)"-^".'- 

(16) 

Here z is the redshifted Lya wavelength (Ac z — 1216[1-|- 

.]A). ' 

In Figure |3 we plot tuab as a function of observed 
wavelength, A in the range A = 2000 — 10000 A for three 
values of the total ffux in the Lya line, /a, 17 ~ 14 {black 
solid line), 2 (red dashed line) and 0.4 {blue dotted line), all 
at 2 = 4.5. We assumed the Lya equivalent width to be 
EW= 140 A. The break at A = 6700 A is the Lyman break 
caused by the IGM; the Lya forest reduces the mean flux 
blueward of the Lya line by a factor of (e"'^). This factor 
has been det ermined observa tionally and is {e~^) = 0.25 at 
z = 4.5 fe.g. lFan et"alll2002i) . Figure 01 (and Eq. 18) shows 
that AGN with fa ^/max = 1.4 x 10^^^ ergs s'^ cm~^ 
have a Bw apparent magnitude of rriAB ^26.6 (A ~ 4400 
A). These AGN would have been detected at the J>5 — a 
level in the Bw band of the NOAO deep wide-field survey 
( J annuzi fc Dc v 199d this detection limit is indicated by the 
thick horizo ntal line at mAB = 26 .6). Since one of the selec- 
tion criteria [Pawson et alj (|20^) use to select z = 4.5 Lya 
emitter candidates states that a z = 4.5 candidate must not 
be dete cted in the Bw b and of the NOAO deep wide-field 
survey jMalhotra fc Rhoa ds 2002). AGN with fa ^/max = 
1.4 X 10~^^ ergs s~^ cm~^ would not have made it into the 
sample. Therefore, with the current selection criteria, the 
Lya survey performed by Dawson ot al. (2004) only probes 
AGN in the luminosity range Lb G [I/B,min, iB.max], with 
Lb, max ~ 2 X lO^^Ls,0 (this luminosity was used in 8^ 13.2.21 . 



4.3 X-Ray Emission from Faint AGN 

X- Rays provide a reliable pointer to AGN activity. Using 
the ISazonov et alj (JiooJ) template, we calculate the ratio 
of detectable Lya and X-Ray fiux (in the observed 0.5-10.0 
keV band) 

fx V0.5/Vl00y4/ ^ ' 



8 Mark Dijkstra and J. Stuart B. Wyithe 



25 



D 



30 



f =14 

h-''=2 1 

'a,17 "-^ 1 

j 




\ 


— _ J : 

, - - - 

- - • * 





2000 



3000 



4000 



5000 



6000 



7000 



8000 



9000 



Wavelength A (Angstrom), Observers frcnne. 

Figure 3. Model spectra of faint AGN at 2 = 4.5 with observable Lya flux, /a = 1.4 X 10^"'^'^ {black solid line), 2 X 10""'^'' (red dashed line) 
and 4 X IQ-^* (blue dotted line) ergs s"! cm-^ obtained with Eq.([TC} with EW=100 A. This figure shows that AGN with /„ ^2 X IQ-^'^ 
ergs s~^ cm~^ have a B\y apparent magnitude of rriAS S^26.6. These AGN would have been detected in the NOAO deep field (whose 
detection limit is indicated by the t hick horizontal line) and would therefore have been discarded as being candidate h igh redshift Lyg 
emitters iMalhotra fc Rhoad j|2002D . Furthermore, the AGN with /„ = 4 X 10 ergs s ^ cm which corresponds to lTaniguchi et alj 
j200^'s detection limit, have apparent magnitudes mAB ^28, which makes them difficult to detect at other wavelengths. 



This ratio changes by less than ~ 10% between z = 4.5 — 6.5. 
The total X-Ray flux in the observers' 0.5-10.0 keV band is 
~ 20 times higher than the total Lya flux. For example, 
an AGN with a detected Lya flux of 2 x 10"^^ ergs s^^ 
cm^^ is expected to have an X-ray flux of ~ 4 x 10^^" ergs 
s~^ cm~^. This lies slightly below the quoted 3— a detection 
limit in the X-ray observations of individual Lya emitter 
candidates presented by Wane, et a l. ( 2004) (~ 6.6 x 10^^'' 
ergs s^^ cm~^ in the 0.5-10.0 keV band). Therefore, even 
if AGN were in the sample, these X-ray observations would 
not necessarily have revealed them ("S 13.1^ . Also, provided 
the fraction of AGN is /agn <,0-^ I3.3II . these would not 
have appeared in the stacked X-ray image, since stacking 
101 images increases the noise by 101^''^, whereas the X- 
Ray signal from AGN would increase by /agn x 100 ;^10. 
The signal-to-noise ratio would thus be conserved at best. 

Similarly, the total X-Ray flux in the observers' 0.5- 
2.0 keV band is ~ 4 times higher than the total Lya flux. 
It may come as a surprise then, that our constraints are 
fainter than those inferred from the Chandra Deep Field. 
The X-Ray detection threshold in the Chandra Deep Field 
Afori/i in the 0.5-2.0 keV band is ~ 1.5x10"" ergss"^ cm"^, 
which would correspond to a Lya flux of ~ 3 x 10~^* ergs s^^ 
cm~^, well below the sensitivity limit of /um ~ 1.1 x 10~^^ 
ergs s ~^ cm~^ in the survey performed bv iDawson et all 
i2004h . However, estimates of the photometric redshifts for 
the X-Ray sources require that they be bright enough to be 
detected in the Subaru Sloan z' band (which is centered on 
A ~ 9000 A). This imposes a flux limit corresponding to AB 
magnitudes of the X-Ray sources in the Sloan z' band that 
are below ^25.2. Using the quasar template, we find that 
this excludes AGN with B-Band luminosities less than 6 x 



10 Ls,Q, which agrees very well with the lower luminosity 
bound of the 'faint' X-ray bin in 13.2.21 

The previous discussion demonstrates an observational 
bias against identification of high redshift faint AGN in deep 
X-Ray observations. Within narrow redshift windows, this 
bias could be alleviated by combining X-Ray observations 
with deep Lya observations, since the detection of a Lya 
line would determine the AGN's redshift. Similarly, the bias 
against identifying bright AGN in wide field Lya surveys 
(§ 14.21 could be alleviated in combination with deep X-ray 
observations. These points illustrate the utility of combining 
deep X-Ray and Lya observations to identify faint AGN, 
provided the sensitivities in each observation probe the same 
population of AGN. 



5 DISCUSSION 

We have shown that ai z = 4.5 and below Mb ~ —20, 
the quasar B-band luminosity function rises more slowly 
towards lower luminosities, dlog'i' /dlog Lb ^— 1.1 (95% 
confidence level), than has been observed at higher lumi- 
nosities and lower redshifts, where 9 log ^'/c' log Lb ~ — 1.6 
(Pci 1995; Bovlc et al. 2000). This fiattening of the faint 
end of the luminosity function to wards highe r reds hift is 
consistent with the recent work bv iHunt et alJ (|2003), who 
found that 91og */91og Li3= -1.24±0.07 at z = 3. We find 
marginal evidence (75% confidence level) that, in fact, the 
luminosity function falls towards lower luminosities below 
Mb ~ —20. This finding is in contrast to observations at 
higher luminosities. Moreover the observed number counts 
lie well below the model predictions of Wyithe & Loeb 2003 
(Fig. 0. This may be explained in three ways: 
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1) Our work has focused on the luminosity function of 
broad-hned AGN (Type I), which in the unified model for 
AGN are the same as Type II AGN (e.g. Norm an et al, 
2002), but unobscured by the thick absorbing torus. IsimpsorJ 
l|2005|) has shown that the fraction of type I AGN increases 
with luminosity, which is supported theo retically by the 're- 
ceding torus' model fe.g. lLawrencellT99ll) . This would imply 
that a luminosity function that incorporates both type I 
and II AGN does not flatten as much at z = 4.5 as shown 
in Figure However, to fully explain the observed flatten- 
ing of the luminosity function requires the fraction of type 
1 AGN to decrease by ~ 1 — 2 orders of magnitude between 
log[Ls/Ls,0] = 10 and ~ 11.5 . Since this luminosity de- 
pendence of the type I AGN fraction is much stronger than 
has been observed, this is very unlikely. 

2) Gas accretion onto black holes in the mass-range 
A/bh = 10^ — lO'^M© is suppressed. One origin of this 
suppression may be negative AGN feedback. AGN with 
log[I/B/I/B,©] = 10.5 are typically powered by black holes 
in the mass range Mbh ~ 10® Mq to a few IO^Mq (§|13J- 
According to the relation between A^bh and the circular ve- 
locity Hcirc of the dar k matt er halo that hosts the black hole 
[Ferrarese fc Merritd J2000l) : iFerraresd (|^02)], this corre- 
sponds to Vcirc = 100 — 250 km s~^. For comparison, suppres- 
sion of accretion due to a photoionised IGM at z = 4.5 — 6.5 
only occurs at fcirc = 40 — 60 km s~^ ('e.g. iDiikstra et'al] 
l2004l) . However. lOekel fc Silk! dlQSd) have shown that super- 
nova driven gas loss as a result of the first burst of star for- 
mation becomes significant in halos with Vdrc ^100 km s^^. 
The latter feedback mechanism would therefore provide a 
more plausible explanation for the reduced gas accretion ef- 
ficiency onto black holes in the mass range Mbh ~ lO^Mo- 
a few W'^Mq. 

3) The number of black holes in the range Mbh = 
10^ — 10^ A/0 is lower than expected f rom the Mbh — Vdrc - 
relation by ~ two orders of magnitude. iHaiman et aP (l2004l) 
used the rareness of black holes with masses Mbh ;^1O''M0 
as a possible explanation for their model of the luminosity 
function of radioloud quasars to overpredict the abundance 
of faint radio sources, by one-two orders of magnitude. The 
possible reduction in the number of black holes with masses 
Mbh ;^1O^M0, may reflect the existence of a minimum su- 
per massive bl ack hole mass, as envis ioned in some formation 
scenarios (e.g. iHaehnelt et al.lll998l . in which the minimum 
black hole mass would be Mbh = W^Mq). 

Our finding that a,t z — 4.5 the QLF fiattens more 
than previously believed for log[LB / Lb,q] ^11 also im- 
plies that faint AGN contribute less photons to the ionising 
background than previously thought. However, this is only 
a small effect since the total ionising photon output from 
AGN per unit volume is oc 

For f3i < 2, this integral is dominated by luminous AGN. 
A more intriguing implication concerns miniquasars, which 
are quasars powered by black holes in the mass range 
Mbh = 10— 10^ M0. It has been suggested that miniquasars 
could contrib ute significantly to the ionising background a t 
high redshift ('M adau et alll2004l : iRicotti fc Ostrikerll2004l) . 
However, if quasar activity decreases with decreasing black 
hole mass (as our results suggest), and if this trend con- 
tinues into the miniquasar-realm, then it follows that mini- 
quasars would not be efficient producers of ionizing radiation 
and would not have contributed significantly to the ionizing 



background. It is worth emphasising that this miniquasar- 
realm may be accessible with existing deep spectroscopi c 
surveys of gravitationally lensed regions iSantos et alj2004) . 
In these surveys, black holes with masses of Mbh ~ 10* Mbh 
can be detected, provided these are accreting at their Ed- 
dington limit fij |4.H . Currently, the only observational con- 
straints on the abundance of high redshift miniquasars are 
derived from the cosmic X-Ray and infrared backgrounds 
(Dijkstra et al. 2004b, Salvaterra et al. 2005). 



6 CONCLUSIONS 

Recent Lya surveys have detected Lya emitting objects 
from redshifts as high as z = 6.5, and at luminosities as 
low as 10 erg s" fe.g. lSantos et alJl2004^ . No evidence of 
AGN activity exists among these several hundred Lya emit- 
ters jPawson et al ] |2004l : IWang et aljl2004l : iTaniguchi et alJ 
[200^. Wide field Lya surveys are designed to deeply image 
wide fields on the sky, yielding survey volumes in a narrow 
shell of redshift space as large as 10^ - 10*^ Mpc^, while 
deep spectroscopic surveys of gravitationally lensed regions 
probe deeper into smaller volumes. The absence of AGN 
within these fields can place a tight upper limit on the num- 
ber density of AGN with Lya luminosities exceeding the 
surveys detection thresholds, I/Lya.c. 

In § 121 we have shown empirically that the Lya lumi- 
nosities of AGN equal their B-band luminosities to within a 
factor of a few (Eq.|lJ. As a result, deep Lya surveys can be 
used to obtain upper limits on the number density of AGN 
with B-Band luminosities exceeding I/B,min ~ 1.4LLyQ,c. 
When expressed in B-band solar luminosities, I/B,min is 
~ 1O**'^Lb,0 (which corresponds to an absolute magnitude 
of Mb = —16). In g 14.11 we demonstrated that such AGN 
are expected to be powered by black holes with masses in 
the range A/bh = 10'' - lO'^M© at z = 4.5 - 6.5. 

We derive upper limits on AGN number densities 
and constrain the quasar B-band luminosity function 
^'(/^B.min, 2) at a luminosity Z/B,min. The non-detection of 
AGN among z = 4.5 — 6.5 LAEs rules out the model predic- 
tions by Wyithe & Loeb (2003), which succesfuUy reproduce 
the brighter end of the observed quasar luminosity function 
at 2: = 2 — 6, at ^ 95% confidence levels at all redshifts. At 
z = 4.5, we find that 91og^'/91ogL_g > —1.6, the value ob- 
served at lower redshifts, for log[L_B/Ls,0] ^11 at the 98% 
confidence level (Figl^J. We find marginal evidence that at 
these luminosities, the luminosity function rises with lumi- 
nosity, corresponding to a powerlaw slope > (75% con- 
fidence level). In other words, the QLF may increase with 
Lb at these faint luminosities, in contrast to observations 
of more luminous AGN. These results represent the faintest 
observational constraints on the quasar luminosity function 
at these redshifts to date. 

We have found that models of the quasar luminosity 
function which are successful in reproducing the bright end 
of the quasar luminosity function predict more AGN to be 
present than are observed, by up to two orders of magnitude 
(Fig. at 2 ~ 4.5. These results imply either that accre- 
tion onto lower mass black holes is less efficient than onto 
their more massive counterparts, or that the number of black 
holes powering AGN with Mb ^ — 20 is lower than expected 
from the Mbh — (y relation by one-two orders of magnitude. 
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Extrapolating from reverberation-mapping studies suggests 
that these black holes would have Mbh = 10*^ - lO'^M©. 

Our work has demonstrated the effectiveness of Lya 
surveys in constraining the faint end of the quasar B-band 
luminosity function. Deeper and larger surveys will allow 
for a better determination of its slope, and whether in- 
deed the quasar luminosity function rises with luminosity 
for Mb ^ — 20 at z = 4.5, and at other redshifts. These 
constraints will offer new insights on the growth of low mass 
black holes and their relation to the known super massive 
black holes. To help identify AGN among observed Lya 
emitters, we have modeled the observable properties of the 
Lya line for high redshift, faint AGN. Using the empirical 
Kaspi-relations, we estimate that the observable Lya line 
widths (Half Width at Half Maximum) of faint AGN will be 
~ 1500 km s"^ For AGN embedded in a neutral medium the 
peak of the Lya line is redshifted up to 2000 km s~^ relative 
to the true line center (ii ld.lll . To facilitate the identification 
of these faint AGN, we have estimated their observable prop- 
erties in the visible (§|lill FigOJ and X-Ray bands (glOl. 
We caution that selection criteria used in Lya surveys to 
select candidate high redshift Lya emitters currently intro- 
duce a bias against detecting AGN with \og[LB /Lb,q] ^11-3 
(i; 14.21 . corresponding to the faintest AGN identified in the 
Chandra Deep Fields. 
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APPENDIX A: TRANSMISSION OF LYa 
PHOTONS FROM AGN 

For a source of UV radiation embedded in a neutral IGM the 
emission blueward of the Lya line center is suppressed by 
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a factor of (e ^'^^). Here, the Gunn-Peterson optical depth, 
TCP is given by 



TGP 



3n{z)AaXi 
8-kH{z) 



= 6.4x10" 



fQth'^ \ ( 0.15 



V 0.022 / Vfi^/i^y V 7.5 

(Al) 

where — 6.25 x 10* s ^ is the Einstein A coefficient, = 
1216 A is the Lya wavelength, and n{z) and Hiz) are the 
number density of neutral hydroge n atoms and the Hubbl e 
constant at redshift respectively l|Gunn fc Petersonll96!Tll . 
For photons initially redward of the Lya line center, the 
Gunn-Peterson optical depth reduces to 



TGP (2;) = TGP —;= 



-TGP- 



10 ( ) V 300K j 1/2 / ~38 



7.5 



(A2) 



Here we have expressed the frequency v in terms of x = 
{v-v^)l I^vd-, where Ar/^ = VtuVo/c, and «th = ^/2k^Tjrn^ 
is the thermal velocity of the hydrogen atoms in the gas, 
ks is the Boltzmann constant, T the gas temperature, nip 
the proton mass and vq = 2.47 x 10^^ Hz is the central 
Lya frequency. Note that a; < for photons redward of the 
line center. To obtain the simple expression in Eq. (IA2II . 
we approximated the Voigt function <j>{x) in the line wing 
as (/"(x) = a/[y/Tvx'^\, where a is the Voigt parameter and 
a = Aa/iirAuo = 4.7 x 10"'' (13 km s"Vuth) is the ratio 
of the Doppler to natural line width. 

Another way to write the Gunn Peterson damping wing 
optical depth is in terms of a velocity offset, Av. For a pho- 
ton initially redshifted by Av relative to the line center, the 
total Gunn-Peterson optical depth reduces to: 



10 



1 -f 2\3/2 /85 km s~i 



7.5 



Av 



(A3) 



Note that this expression is independent of the gas temper- 
ature. 

Assuming the intrinsic Lya spectrum of an AGN to be 
Gaussian of width CTv = hhwhm (Eq. 1131 with a central ffux 
that is A'^ times the continuum''. The intrinsic flux density 
(in arbitrary units) becomes 



N-l r-X^. 

l + ^^exp[^]. 



where = o^/vth- The IGM transmission is given by 



(A4) 



/_°°^dx[F(x)-lle- 
dx[F{x) - 1] 



-r(x) 



(A5) 



where F{x) and t{x) are given by Eq. l|X4j and Eq. llX2t . 

The transmission T is plotted as a function of in ¥ig- 
ure lAll as the hlack solid line. 

The other lines represent the IGM transmission in the 
presence of a cosmological HH region with a radius of 0.1 
[red-dotted line), 0.5 [green-dashed line) and 1.0 Mpc [blue- 
long-dashed dotted line) surrounding the AGN. The radius 
of the Hn region surrounding a quasar of age tq powered by 



^ This corresponds to a Lya EW of N X [Sijhwhm/c] x Aq 
A=200 X [N/10) X (t)HWHM/[3000 km s^i]) A. 



R„,i-1 Mpc 
R„„-500 kpc 
R„,i 100 kpc 
No HJl region 



Figure Al. The fraction T of Lya photons emitted hy a z = 6.5 
quasar embedded in a neutral IGM, that is transmitted to the ob- 
server, as a function of the intrinsic Lya line width of the quasar, 
(Tv = whwhm- The four lines represent models with cosmologi- 
cal HII regions of different sizes. The figure shows that an IGM 
transmission of T = 0.3 is accurate to within a factor of 2 for a 
wide range of parameters. 



a black hol e of mass Mbh shining at Eddington luminosity 
IS given by JHaiman fc Cen,2002i 



Rmi.t 



= 0.74 



Mb 



IO^Mq 



1/3 



2 X lO^yr 



1/3/1+2 



7.5 



Mpc. 

(A6) 

This relation assumes the total ionizing luminosity of AGN, 
Lion, to scale as Lion oc M. Alternatively, the Kaspi relation 
shows that the continuum luminosity at A = 5100 A scales 
as oc M^ '^ (Eq. [TTJ- H the brightest z > 6 AGN accrete 
at the Eddington luminosity, and we assume that the total 
ionizing luminosity emitted by an AGN also scales as M^'"^, 
we get 



RniiMs = 0.44 



Mbh 
lO^Mo 



_tQ y/^n+z 



2 X 



W^yrJ 



7.5 



Mpc. 

(A7) 

The size of the HII region is also determined by the 
uncertain value of the escape fraction of ionizing photons, 
/esc. Accounting for these uncertainties. Figure IXTl shows 
that for av — 3000 km ,T lies in the range 0.3—0.45 given 
Rhii = — 500 kpc. The range of T increases towards lower 
(T^. However we find that an IGM transmission of T = 0.3 
is accurate to within a factor of 2 for a wide range of the 
parameters a^, /cse and total ionizing luminosity. 

In Figure IX2I the black solid line shows the theoretical 
spectrum of a z = 6.5 AGN embedded in a neutral IGM (i.e. 
no cosmological HII region is present). We assume the black 
hole mass to be lO' Mq, which according to Eq. 11311 should 
yield whwhm = 4x 10'' km s"'^. However, at low luminosities, 
and therefore low blac k hole masses, whwhm — > 3 X 10^ km 

iKaspi et alJl200oli . The intrinsic Lya spectrum (the red 
dotted line) is described by Eq. IIA4L in which we assumed 
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Faint AGN at z=6,5 (x„,"1.0) 




8800 9000 9200 9400 9600 9800 

A(A) 



Figure A2. The observed spectrum of a faint AGN at 2 = 6.5 
in a neutral IGM (i.e. no cosmological HII region is present) is 
shown as the black solid line. The intrinsic Lyo line is shown as 
the red dotted line. A neutral IGM erases a large fraction of the 
line, transmitting only the reddest photons to an observer. The 
figure shows that these AGN have a peak Lya flux redward of 
the true line center by ~ 60 A(corresponding to ~ 2000 km s~^) 
and an observed f hwhm that is a factor of ~ 2 times lower than 
the intrinsic fHWHM. 



o-„ = 3 X 10^ km s"^ and iV = 10 (which yields an Lya EW 
of 200 A). 

We note that the observed line center lies redward of the 
true line center by ~ 60 A, which translates to ~ 2000 km 
s~^. Furthermore, the observed «hwhm is 2 times lower 
than the emitted value. 



